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Abstract 



C \ In this letter, taking the weU known (2+l)-diniensional soHton systems, Davey-Stewartson 

(DS) model and the asymmetric Nizhnik-Novikov-Veselov (ANNV) model, as two special 
examples, we show that some types of lower dimensional chaotic behaviors may be found in 
\^ , higher dimensional soliton systems. Especially, we derive the famous Lorenz system and its 

0^ . general form from the DS equation and the ANNV equation. Some types of chaotic soliton 

(N ■ 

f^ . solutions can be obtained from analytic expression of higher dimensional soliton systems and 

r_^ . the numeric results of lower dimensional chaos systems. On the other hand, by means of 

the Lax pairs of some soliton systems, a lower dimensional chaos system may have some 
types of higher dimensional Lax pairs. An explicit (2+l)-dimensional Lax pair for a (1+1)- 
dimensional chaotic equation is given. 

^. ' 

.^ ■ PACS numbers: 02.30.Ik, 05.45.-a, 05.45. Ac, 05.45. Jn 

x: 

H . In the past three decades, both the solitons |]l| and the chaos||2[ have been widely studied and 

applied in many natural sciences and especially in almost all the physics branches such as the 
condense matter physics, field theory, fluid dynamics, plasma physics and optics etc. Usually, 
one considers that the solitons are the basic excitations of the integrable models, and the chaos 
is the basic behavior of the nonintegrable models. Actually, the above consideration may not 
be complete especially in higher dimensions. When one says a model is integrable, one should 
emphasize two important facts. The flrst one is that we should point out the model is integrable 
under what special meaning(s). For instance, we say a model is Painleve integrable if the model 
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possesses the Painleve property, and a model is Lax or 1ST (inverse scattering transformation) 
integrable if the model has a Lax pair and then can be solved by the 1ST approach. An integrable 
model under some special meanings may not be integrable under other meanings. For instance, 
some Lax integrable models may not be Painleve integrable |^, ^. The second fact is that for the 
general solution of a higher dimensional integrable model, say, a Painleve integrable model, there 
exist some lower dimensional arbitrary functions, which means any lower dimensional chaotic 
solutions can be used to construct exact solutions of higher dimensional integrable models. 

In this letter, we will show that an 1ST integrable model and/or a Painleve integrable model 
may have some lower dimensional reductions with chaotic behaviors and then a lower dimensional 
chaos system may have some higher dimensional Lax pairs. 

To show our conclusions, we use the (2+l)-dimensional Davey-Stewartson (DS) equation|^] 



iut + 2 {uxx + Uyy) + a\u\ u — uv = 0, (1) 



Vyy - 2a{\u\'^)xx = 0, (2) 



as a concrete example at first. The DS equation is an isotropic Lax integrable extension of the 
well known (l+l)-dimensional nonlinear Schrodinger (NLS) equation. The DS system is the 
shallow water limit of the Benney-Roskes equation^], where u is the amplitude of a surface 
wave-packet and v characterizes the mean motion generated by this surface wave. The DS 
system (1) and (2) can also be derived from the plasma physics|^] and the self-dual Yang-Mills 
field. The DS system has also been proposed as a 2+1 dimensional model for quantum field 
theory]^]. It is known that the DS equation is integrable under some special meanings, namely, it 
is 1ST integrable and Painleve integrable Q. Many other interesting properties of the model like 
a special bilinear form, the Darboux transformation, finite dimensional integrable reductions, 
infinitely many symmetries and the rich soliton structures [P, ^ have also been revealed. 

To select out some chaotic behaviors of the DS equation, we make the following transforma- 
tion 

V = Vo- rHfx'x' + fy'y' + 2/xv) + f'Hf!' + Vy' fx' + §), (3) 

u = gf-^ + uo (4) 

with real / and complex g, where x' = {x + y)/\/2, y' = (x — y)/\/2, and {uq, vq} is an arbitrary 
seed solution of the DS equation. Under the transformation (3) and (4), the DS system (1) and 
(2) is transformed to a general bilinear form: 

{Dx'x' + Dyiyi + 2iDt)g ■ f + uq{Dx'x' + 2Dx'y' + Dy/y/)f ■ f 

+2auogg* + 2a4g* f - 2vogf + Gifg = (5) 



2{D,,y, + a|noP)/ • / + 2agh + 2agful + 2auog* f - GJf = 0, (6) 



where D is the usual bilinear operator|Tn] defined as D^A ■ B = [dx — dxi)^A{x)B{xi)\x'i^=x, 
and Gi is an arbitrary solution of — 16a(tiox' + '^Oy')Wox' + ^Ov') + Gix'x' + Giy'y' + 2Gix'y' — 
4:a{Dx'x' + Dyiyi + 2Dx'y')uQ ■ Uq = 0. For the notation simplicity, we will drop the "primes" of 
the space variables later. 

To discuss further, we fix the seed solution {uq, vq} and Gi as 

UQ = Gi = 0, vo=po{x,t) +qo{y,t), (7) 

where po = po{x,t) and qo = qo{y,t) are some functions of the indicated variables. 
To solve the bilinear equations (5) and (6) with (7) we make the ansatz 

f = G + p + q, g = piqiex.p{ir + is), (8) 

where p = p{x, t), q = q{y, t), pi = pi{x, t), qi = qi{y, t), r = r(x, t), s = s{y, t) are all 
real functions of the indicated variables and C is an arbitrary constant. Substituting (8) into 
(5) and (6), and separating the real and imaginary parts of the resulting equation, we have 

2pxqy - apjql = 0. (9) 

(qiPixx + PiQiyy - Piqii^n + 2st + 2{po + go) + 4 + ''x))(C + p + q) 
+qi{piPxx - 2pixPx) + piiqiqyy - "^qiyqy) = (10) 

i-qi{2rxPix + 2pit+pirxx) - Pi{2syqiy + 2qit + qiSyy)){G + p + q) 

+2qipi{qt + Syqy) + 2qipi{rxPx + Pt) = (11) 

Because the functions p^, p, pi and r are only functions of {x, t} and the functions qo, q, qi 
and s are only functions of {y, t}, the equation system (9), (10) and (11) can be solved by the 
following variable separated equations: 



pi = 5iy2a-iq V', qi = hy/ci%, {5l = 6l = 1), (12) 

Pt = -rxPx + C2, qt = -Syqy - C2, (13) 

4{2rt + rl + 2po)pl + pL - '^PxxxPx + cqpI = 0, (14) 

4(2sj + 4 + '^qo)ql + qyy - "^qyqyyy - coqy = 0. (15) 



In Eqs. (12)-(15), Ci, C2 and cq are all arbitrary functions of t. 

Generally, for a given pQ and qq the equation systems {(12), (14)} and {(13), (15)} may not 
be integrable. However, because of the arbitrariness of po and qo, we may treat the functions 
p and q are arbitrary while po and qo are determined by (14) and (15). Because p and q are 
arbitrary functions, in addition to the stable soliton selections, there may be various chaotic 
selections. For instance, if we select p and q are solutions of (ri = x + uJit, T2 = x + L02t) 

PriTin = [PtitiPti + (c + l)p?J P'^ - {p"^ + bc + b)pri - {b + c + l)priri + pc{ba -b- p^), (16) 
<1t2T2T2 = [qT2T2QT2 + (7 + 1)9?2) 1^^ " (^^ + Pi + P)<1t2 " (/^ + 7 + l)gr2T2 + qc(l3a - (3 - 9^), (17) 

where uji, 102, Oj b, c, a, /3 and 7 are all arbitrary constants, then 

Co = C2 = 0, r = -L(Ji(x + u;it/2), s = -u)2{y + "^2^/2), (18) 

Po = -4"^ (cpV/ + 15^ - bcp-^ia - l)p + b{c + l) + {b + c+ 1)^^,^,^-^ 

+pl^^^2-^p-^-p-\p,^{c+l)+p^,r,)) , (19) 

90 = -4-^ (79^97/ + 9^ - /?79r2"^(a - 1)9 + /3(7 + 1) + (/? + 7 + I)gr2r2 9r/ 

+d 2-ig~2 _ q-\q^^[^ + 1) + g^^^j) . (20) 



Substituting (8) with (12)-(20) into (3) and (4), we get a general solution of the DS equation 



u = 5i52\j2a iprigr2exp(-i(a;i3; + a;2y+ 2(^1 +^2)*)(C + P + Q') ^ (21) 

v=Po + qQ- {qr2T2 + Pr,T,){C + p + q)-^ + {ql^ + 2qr^Pr, + pI,){C + p + q)-^ (22) 

where pq and go are determined by (19) and (20), while p and q are given by (16) and (17). 

It is straightforward to prove that (16) (and (17)) is equivalent to the well known chaos 
system, the Lorenz system[ll|: 



Pti = -c{p - g), 5ri = Pia - h) - g, h^^ = pg - bh. (23) 

Actually, after canceling the functions g and h in (23), one can find (16) immediately. 
Prom the above discussions, some interesting things are worth emphasizing: 
Firstly, because of the arbitrariness of the functions p and g, any types of other lower dimen- 
sional systems may be used to construct the exact solutions of the DS system. 

Secondly, the lower dimensional chaotic behaviors may be found in many other higher di- 
mensional soliton systems. For instance, by means of the direct substitution or the similar 
discussions as for the DS equation, one can find that {ji = x + uoit) 

u = 2p^Wy{p + w)-'^, (24) 



2Px , ic + l)Px {bp - w)p,j,,j, 1 2 1 rn , ^ , ^ , 3 



■^ , ^ ,2 + ^ - o , ' \ - o^ - :^[(1 + fc + cW + cp' - cb{a - l)p], (25) 

[p + w)'^ dp 6p[p + w) 6 6px 

with w being an arbitrary function of y, and p being determined by the (l+l)-diniensional 
extension of the Lorenz system 

Pt = -Pxxx + P~^ \PxxPx + (c + 1)pI] - p^Px - (6 + c + l)pxx - pc{h -ba+p^), (26) 

solves the fohowing 1ST and Painleve integrable KdV equation which is known as the ANNV 
model m 

ut + Uxxx - Huv)x = 0, Ux = Vy. (27) 

It is clear that the Lorenz system (16) is just a special reduction of (27) withp = p{x+b{c+l)t) = 
p{ti). Actually, p of equations (24) may also be arbitrary function of {x, t} after changing (25) 



appropriately |12]. In other words, any lower dimensional chaotic behavior can also be used to 
construct exact solutions of the ANNV system. 

The third thing is more interesting. The Lax pair plays a very important and useful role in 
integrable models. Nevertheless, there is little progress in the study of the possible Lax pairs 
for chaos systems like the Lorenz system. In Ref. [Q, Chandre and Eilbeck had given out the 
Lax pairs for two discrete non-Painleve integrable models. In Ref. 0, we have found some Lax 
pairs for some non-integrable Schwarzian equations. Now from the above discussions, we know 
that both the Lax pairs of the DS equation and those of the ANNV system may be used as 
the special higher dimensional Lax pairs of arbitrary chaos systems like the Lorenz system (16) 
and/or the generalized Lorenz system (26) by selecting the fields appropriately like (21)-(22) 
and/or (24)-(25). For instance, the (l+l)-dimensional generalized Lorenz system (26) has the 
following Lax pair 

ipxy = uip, ipt = -tpxxx + 3vipx (28) 

with {u, v} being given by (24) and (25). From (24)-(27), we know that a lower dimensional 
chaos system can be considered as a consistent condition of a higher dimensional linear system. 
For example, ipxyt = i^txy of (28) just gives out the generalized Lorenz system (26). 

Now a very important question is what the effects of the lower dimensional chaos to the 
higher dimensional soliton systems are. To answer this question, we use the numerical solutions 
of the Lorenz system to see the behaviors of the corresponding solution (24) of the ANNV 
equation by taking w = 200 + tanh(y — yo) = Ws and p = p{ti) = p{X) as the numerical solution 
of the Lorenz system (16). Under the selection w = Ws, (24) is a line soliton solution located 
at y = yQ. Due to the entrance of the function p, the structures of the line soliton become 




Figure 1: Plot of the period two line soliton solution (24) of the ANNV equation with (29), 
P = pi''')-, aiid a = 350, b = 8/3, c = 10 at t = 0. 



very complicated. For some types of the parameters, the solutions of the Lorenz system have 
some kinds of periodic behavior, then the line soliton solution (24) with w = Wg becomes an x 
periodic line soliton solution that means the solution is localized in y direction (not identically 
equal to zero only near y = yo) and periodic in x direction. Fig. 1 shows the behavior of the 
periodic two line soliton solution and p being a periodic two solution of the Lorenz system (16) 
with a = 350, b = 8/3 and c = 10. In some other types of the parameter ranges, the solution 
of the Lorenz system becomes chaotic and then the line soliton of the ANNV system becomes 
chaotic line soliton which is localized in y direction and chaotic in x direction. Fig. 2 displays 
the chaotic behavior of the amplitude of the line soliton located at y = yo with a = 60, b = 8/3 



and c = 10. The parameters we used here are the same as those used in literature 1 11]. 

In summary, though some (2+l)-dimensional soliton systems, like the DS equation and the 
ANNV equation, are Lax and 1ST integrable, and some special types of soliton solutions can 
be found by 1ST and other interesting approaches H, any types of chaotic behaviors may still 
be allowed in some special ways. Especially, the famous chaotic Lorenz system and its (1+1)- 
dimensional generalization are derived from the DS equation and the ANNV equation. Using 
the numerical results of the lower dimensional chaotic systems, we may obtain some types of 
nonlinear excitations like the periodic and chaotic line solitons for higher dimensional soliton 
systems. On the other hand, the lower dimensional chaos systems like the generalized Lorenz 
system may have some particular Lax pairs in higher dimensions. 

It is also known that both the ANNV system and the DS systems are related to the 
Kadomtsev-Petviashvili (KP) equation while the DS and the KP equation are the reductions 
of the self-dual Yang-Mills (SDYM) equation. So both the KP and the SDYM equations may 




Figure 2: Plot of the chaotic hne sohton solution (24) of the ANNV equation with the same 
conditions as Fig. 1 except for a = 60. 

possess arbitrary lower dimensional chaotic behaviors induced by arbitrary functions. From the 
results of this paper, various important and interesting problems should be studied further. For 
instance, what on earth is the complete integrability and what kinds of information about chaos 
can be obtained from some types of special higher dimensional Lax pairs? 
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